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Abstract 

Dixit and Jabbari Nooghabi (2010) had derived the MLE and UMVUE of the probability density 
function (pdf) and cumulative distributive function (cdf). Further, it had been shown that MLE is 
more efhcient than UMVUE. He, Zhou and Zhang (2014) have also derived the same and made a 
remark that the work of Dixit and Jabbari Nooghabi (2010) is not correct. We have made a comments 
with detail algebra that our results are correct. Further, we have also given the R code. 
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1 Introduction 

The Pareto distribution has been used in connection with studies of income, property values, insurance 
risk, migration, size of cities and firms, word frequencies, business mortality, service time in queuing 
systems, etc. 

The objective of this paper is to discuss efficient estimation of pdf and CDF of Pareto distribution which 
has been one of the most distinguished candidates for the honor of explaining the distribution of incomes, 
assets, etc. 

We assume that random variable X has Pareto distribution with parameter a and k (known) and its 
probability density function (pdf) is as, 

Oik^ 

fx{x) = 0<fc<a;, a>0 

and distribution function ^ 

Fx{x) = l-(^^^ , k<x 

In economics, where this distribution is used as an income distribution, k is some minimum income with 
a known value. Asrabadi (1990) derived the uniformly minimum variance unbiased estimator (UMVUE) 
of the probability density function (pdf), the distribution function (cdf) and the moment. 

In this paper, we will give the detail algebra of Dixit and Jabbari Nooghabi’s (2010) paper. Also, We 
have made a comments that our results are correct. Further, we have also given the R code. 

Dixit and Jabbari Nooghabi (2010) had derived the MLE and UMVUE of the probability density function 
(pdf) and cumulative distributive function (cdf). Further, it had been shown that MLE is more efficient 
than UMVUE. He, Zhou and Zhang (2014) have also derived the same and made a remark that the work 
of Dixit and Jabbari Nooghabi (2010) is not correct. 

We like to make some comments as follows. 

1. We have verified our results and they are correct. We have given the detail algebra and R program. 
See the attachment. 

2. Examples given by He et al. (2014) are not correct. One should note MSE of f{x) or F{x) is a function 
of parameters. By this one cannot prove anything. Only, one can calculate /(x) or F(x). 
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3. According to definition of the modified Bessel function in Olver, Lozier, Boisvert, et al. (2010) in the 
Theorem 1. the notation K/^n-r)nraz) should be K(^r-n){‘^\/nraz). Also, the notation K(^n){2y/nraz) 
should be K(^_n) {2y/nraz). 

4. After the Theorem 2. the Kummer confluent hypergeometric function is wrong and the correct version 
is 


U (a, &, c) 


1 
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2 Main Result 


In this section, we give the detail algebra of the paper Dixit and Jabbari Nooghabi (2010). 
The details of finding result of the second chapter of that paper are as follows. 


1 ,... 

fix) 

ak°‘ 

- 2 , 0+1 ’ 

L{xi 

, ... , Xyi 1 ct, k) 

a”fc' 

T— 


0 < k < X, a>0 


■ ]^I(x* - k), 


where I is the indicator function defined as 



y>o, 

otherwise. 


( 1 ) 


Put: 


In T(x, a,k) = n ln(a) + na In(fc) — (a + 1) ln(^^ Xi) 


i91nL(s,a) n , ,,, , 

=-l-nln(fc) —ln> Xi = 0 

rv < ^ 


da a 


d = MLE{a) = 


i=l 

d'^lnL -n 


Er=iln(f)’ a 2<0 

dfc" 

O’<5+1 


^ MLE of f{x) = fix,d) = 

~ cxk^ 

/(a;) = 5 > 0, 0<k<x, 


Fix) = 1— ( — I , a>0, 0<A:<x. 

\ X / 


fy(y) = key-^ 

{key) 


- = ae~°‘y, y > 0, x > k 


fY{y) = ae°‘y, y > 0 or r~r(l, i) 


n ^ 

5 = ^y,~r(n,-), 

r\ 


9s{s) = 




exp(—as), s > 0. 


( 2 ) 


( 3 ) 
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Let w = a 


Jacobian = I—^1 = — k- 


,,, n n ds —n 

il)^w = -^s= — ^ — — 

s w aw 

>g{w) = w>0, a>0 

(an)" r an'i 

^ ^ ^ = rt \f n+n ®^P 1-P w > 0. 

i {n)[w^^^) I w ) 


E{a) = E{W) = / wg{w)dw = / w exp -Uw = —— 

Jo Jo r(n)(n;"+^) t w ) r(n) J^ 


w "e ™ fiw. 


Put z = — ^ dz = — ^dw, so 

•3/1 7/1-^ ' 


P(tP) = 


z"-2e-“"^dz = 


-r(n-i)(—)"-i / 

Jo 


oo ^n— 2 g—Q;n 2 ; 


L;(a) = 


^(a^) = = 

Jo 


/•- 2 («n)" f (an)" r ^ 

/o ^ r(n)(n;"+i)®^Pl n; / " r(n) X 


Same as the previous 


= w y„ ' ' ‘'' = -rw"‘" ■"*'»> y, r(„-2)(i)- 


P(a2) = 


(n — l)(n — 2) 


Therefore 


MSE(W) = V(W) + (pytp) - a)^ = pytp^) - 2aE(W) + a^ 

an 2 

/ r)\ -1 E CX . 

(n — l)(n — 2) n — 1 


MSE{a) = MSE{W) = 


+ n — 2) 
(n — l)2(n — 2) 


Proof of the Theorem 1.: 

(A) 


2 2 f ff \ I \.i f (an)" _„„/ 

E(f(x)) = / f{x)g(w)dw = / -n—-rre ' 

y-iy // j jy jyy j ^n,+i r(n)n;"+i 

(an)" /•“ /A:\“ e-“"/“ , 


r(n)x Jo \a ^7 


Put (^)"' = then 


E{f{x)) = 


(an)" r“ 

r(n)x Jo 


3 



We know ™ ^ 


Eim) = 


(an)” f 
r(n)a; Jq 


E 




j=0 

CJO /!„ / k\\j 


jlyjn 


^/^dw 


(an)” ^ (ln(f)) 


^—OLriji} 


r(n)x ^ j! 7o 


-dw. 


Put ^ = z then -^dw = dz ^ dw = -^. 
E{f{x)) : 


(B) 


^n-j-2^-anz 


(an)»^ (ln(|))^ r(n-j-l) _ 

r(n)x ^ j! (^)-"+J+i io r(n -j- l)(;^)”-^-i 

1 ^ (na)J+i 

fW,Z^^nn-,-r,[,n[- 

n—2 


dz 


Eifix)) = 




r(n)x ^ j! 


( 5 ) 


E{F{x)) = / = 


(an)^ 


i-(ir 

X 


{any 


\ r(n)z 


,n+l 




= 1 - 


frvr?'!^ A’ 1 

Q r(n)'u;”+i r(n) Jg w”+i 

(an)” /•“ e”'‘”(-)e-“”/’' 


We know that (f)^ = ^ ^ then 


r(n) Jo w 

(ln(|)) 


n+1 


-dw. 


i=0 j! 

r(n) 7 o ^ j! 

= 1 _ g MIE /■“ n;-”-i+^e-“”/“dn; 


r(n) 


i=o 


J! do 


r(») j'. to »"+■-' 

Put — = 2 zL^dw = dz ^ dw = zz^dz. Then 

00 ^^.^^k\\j / ^ \'>^—j roo ^n—j — l^—anz 


E{F{x)) = l-^f;M^r(n-j) 

1 in) 7! 

^ ' j=o •' 


z"- "e 

any jg r(n - j)(^)”"^' 

3 


dz 


= i-rM£^r(.. -j)(^in(^i)) , i>t. 


( 6 ) 


4 



Proof of the Theorem 2.: 

(A) At first, we should find E{f{x))'^. So 


^2 _ / (an)" 


Eifix)) = / (fix)) g(w)dw = 


(any 


r(n)a 



X / w 


n—1 


>^210-1-2 


-dw. 


r(n)r 


,n+l 


/^dw 


Similarly to the pervious Theorem, we have 
(an)" /•“ 




(ilf; 


It;" 


E 


2W(ln(|))J' _ 




r(n)x^ 

(an)" 

r(?^)a;^ 7o ^ 

(an)" ^ 2J(ln(f))J' ^ 

r(n)x2 j! “ - ^ 

(«n)" - 2^(ln(|))^- 
r(n)x2 ^ -i'l 




z"-^-3e-“"^d0 


i=o 




^ \ /.oo ^n—j —3g—an^; 


r(n)x 


i=o 

n—2 


j! 


r ( n - j - 2)(^)"-^-2 


dz 




r(n) 




x^ 


j=0 


j! 


We know that V(f{x)) = E(f(x)y — E^{f{x)). Then 

n—2 

V{f{x)) = 


Therefore 


(B) Seme as the case (A) 


r(n)x 




j=0 


MSE(f(x)) = 


MSE(f(x)) = 


1 (an)'’+3 . 


V{f{x)) + (Eifix)) - f{x)f 

V(f(x)) + E^fix)) - 2E(f(x))f(x) + f(x) 

E(f(x)f - Eyjix)) + E\f(x)) - 2E(f(x))f{x) + f(x) 
E(f(x)f-2f(x)E(f(x))+f(x) 

1 v 2^(ln(| ))^\ 

r(n)a;2 ^ j! 

nV £ f £ t + 

x^^^l{n)xp^ jl \x J x^^^ 


-r(n — j — 2){any~^‘^ 


E(E(x)f = / {F(x)fg{w)dw = 

Jo Jo 


“ 'l_ 

x'“ 


(any 


T(n)w 


n+1 


"/“dw 
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{any 


T{n)i 


,n+l 


-‘^^l^dw - 2 


(an)" 


roo 

-dw + 


,,n+l 


Tin) 


CO ^ l^\2w „ —an / ^^ 

X ' _ 

.,,n+l 


(an)" /•“ (an)" /•“ 

= '"'wyo 


EiFix))^ = 


1 - 2 

(any 

TR 


(an) 


„ /-oo yj3 


yi-(i)) 


j ^—anjw 


0 


r(n) 

, «00 OO 


/o 


E 

i=0 


E 

j=0 •'' 

2^wyin{^))^ g-an/t» 


■n+l 


dw 


r- 


w 


■n+l 


dw 


i_2feEgEll)l 


j foc ^—anjw 


r(n) ?! 

^ ^ j=o 


/O 


,n+l-? 


(ilf; 


w 


(an)"g2^(ln(f))^ /“e 


oo —OLnjw 


r(nj — J 


,n+l-j 


dw 


w 


“= 1-2 


(«nr^(ln(|)P 


■E 


TR j^o d- 
(«„)n - 2^(ln(|))^- 


r(n-j) 


1 \ ^ ^ /*co ^n—l—jg —0:712 


■E 


r(n) J! 


r(n-i) — 


'o r(n-j)(^) 


2 rco ^n—1—j^ — a.nz 


1 


rdz 


an; do r(n-j)(^) 


1 


rdz 


i?(Rx))2 = 1- 


2 RR(n-j)(an)J 


r(n; 


E 

?=0 


j! 


(In - ) 


1 ^R(n-j)2^ (an)J 


r(n) 


E 

j=o 


(iM-)r 


> k. 


We have V{Fix)) = E{F{x))'^ — E‘^{Fix)), so 


ViFix)) = 1- 


r(n) 


E 

3=0 


r(n — j){an)d 


r- 


(In - y + 


1 ■R r(n — j)2^(an)'’ 


r(n) 


E 

3=0 


r- 


(In(-)F 


1 - 


r(n) 


E 

3=0 


T{n- j)iany 


j! 


(In - ) 


Further MSE{F{x)) = ViF{x)) + iE{F{x)) — F{x))'^ = E{Fix))"^ — 2F{x)EiF{x)) + F'^{x), then 


MSEiFi.)) = 


- 2 


r(n) ^ 

^ ^ 3=0 


1 - 


j! 


1 ^r(n-j)2^(an)^,, 

-7^-(In 


r(n) 


3=0 


j! 


,_r(n-j)(an)^' 

r(n) ^ 

^ ^ 3=0 


f- 


(ln( - ))^ 


1 - 


Therefore 


MSE{F{x)) = 2 + 


r(n) 


E 

3=0 


r(n — j)2^(an)^ 


j! 


(In(-)r 


dw 
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_ 2 I'^V ^ V r(n - J)(»ny (t]y + ^'‘ 

j\ V ^ / V ^ 


X J r(n) 




2oc 


From Asrabadi (1990), we have 


the UMVUE of f{x) and F{x) is 


n — 1 


a = 


ln(t) — nln{k) 


t>k^ 


f(x) = - l)[ln(^) - ln(a:) - (n - 1) In(fc)]" " k<x< tk^-^ 


x[ln(t) — nln(fc)] 


( 8 ) 


(9) 


( 10 ) 


0 


X < k, 


TP<rr.\—} 1 [ln(t)-ln(a;)-(n-l) In(fc)]" ^ l„ ^ ^ ^ 4.lA-n 

X{X) - \ i- [!„(()-„ ln(fc)]'‘-l- ftSX Stti 

1 x> tk^~^. 


( 11 ) 


respectively. Also, 


(xk^ 

fi^) = -^7 x>k>0, a>0, 


F{x) = 1 — J , k < X. 

Proof of the Theorem 3.: 

(A) It is obvious that E{f{x)) = f{x). So, we should find 


E{fix)) = / / ix)h*ix)dt, 


where, 


h*{x) = - - —t “ ^[ln(t) — nln(fc)]” t> fc". 

[n — Ij! 


Therefore 
E{f{x)Y = 


2 _ / (n-l)2[ln(t)-ln(a;)-(n-l)ln(fc)]2" 4 


Ixk^-^ 


a;^[ln(t) — nln(A:)]^” ^ (n—1)! 

(n — 1 )^ 0 ;"fc“” f°° [ln(t) — ln(a;) — (n — 1) ln(fc)]^”“4 


t-“-4[ln(t) -nln(fc)]"-4dt 




x‘^{n — l)l Jxk-^-^ [ln(t) — nln(A:)]” ^ 

Let 2 = ln(t) — nln(A:) ^ dz = \dt, then 

E{f{x)f = ~ n - ln(x) + ln(fc)]^"-4e-^^ 


x^{n-2)\ 


yTl—l 


We know that 


2n-4 


[z - Hx) + ln(fc)]^"-4 = ^ C'(2n - 4, j) (- In (|))' 

i=o 
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where C{n,k) = So 


2n-4 


S C(2n - 4. j) (- l.(f 


fcy 


n-3-j -az 


e-°^^dz. 


3=0 ■' 

The above integral is the incomplete Gamma function, therefore 

2n-4 


^ ^ 3=0 

^ T{n-2-j) "g" exp (-gi n (D) (aln(f))' 


a 


n-2-j 


7.1 


(n — l)a‘^k^ 


2=0 

2n-4 


x^+2(^_2)! 


C{2n-4:,j)a^ 


-In 




2=0 


We know that the Gamma function is defined on the positive value. So 

>' 3=0 

Finally 


2=0 


MSE{f{x))=V{f{x)) = 


(n — l)a^fc° 
a;“+2r(n - 1) 


^C'(2n-4,j)a^r(n-j-2) (-ln(|))^ 

3=0 


"-^-a*(ln(|))^ 


E 

2=0 


rp(y.-\-\ 


(B) 


E{F{x)f = j F'^{x)h*{t)dt 

-h*{t)dt+ / h*{t)dt 


' xk^~ 
roo 


[ln(t) - ln(x) - (n - 1) ln(fc)]2"-2 


pOO 


' xk^ 


[ln(t) - nln(A:)]2"-2 

[ln(t) - ln(x) - (n - 1) ln(fc)]"-i 

k'^ J xk'^~^ 

-a-1 , 


pOC pXK pOO 

/ h*{t)dt+ h*{t)dt-2 

Jxk^~^ J k'^ J xk'^~^ 

t-“-i[ln(t)-nln(fc)]”-idt 


[ln(t) — nln(A:)] 


n—1 


(n — 1)! 


/ xk'^~'^ 

POO 


[ln(t) — ln(a:) — (n — 1) In(fc)]^” ^ 

[ln(t) — nln(A:)]2”“2 (n—1)! 

,,n 7,an poo 


t-“-i[ln(t)-nln(A:)]"-idt 


poo n-LOLYl poo 

/ h*{t)dt — 2- -— / [ln(t) — ln(x) — (n — 1) ln(fc)]"“^t“““^(it 

Jfe" ~ 1)' Jxk"-^ 

[ln(t) — ln(x) — (n — 1) ln(fc)]^"“^ 


a”fc“” f°° 

(n 1)! Jxk^~^ 


[ln(t) — nln(fc)] 


n — 1 


t-°‘-^dt. 


( 12 ) 



We know that h*{t)dt = 1. For second part let 2 ; = ln(t) — ln(a;) — (n — 1) In(fc) and to solve the third 
integral put 2 : = ln(t) — n\a{k). Then 


E{F{x)f = 1-2 


(n — 1)! 




= 1-2 


(n- 1 )! 7in(f) 
a;“(n- 1)! Jg 


[z — ln(a;) + ln(A:)] 

^n —1 


2n —2 


Zu-l^-azdz + 


[z — ln(x) + ln(A:)]^^ 


2n—2^—az 


(n- 1)! 


-dz. 


/p°° z” and for the last integral, we know [z — ln(a:) + In(fc)]^” ^ ~ 

2,j)z^"“^“'l Therefore 


2n—2 




2n—2 




-In 


-In 


(f) 

(!) 


Z«-l-ig-aZzn-lrfz 


i r(n-j) 


^-az 


a”-^' 7i„(|) r(n-j) 


dz. 


The last integral is the incomplete Gamma function, then 

2n-2 


EiFMr = l-2- + ^i:c(2„-2,i)[-ln(|) 


j=0 


rv^ 3 • ^ 7 1 


2=0 


The Gamma function is defined over positive value. So 


n—1 


EiP(z)f = 1-2- + -—i:C(2..-2,i)[-l„(^) 


x^{n — 1)! ^ 

^ ^ i-o 

X a'’r(n-j) ^ 


7 ! 


2=0 


Then 


MSE{F{x)) = V{F{x)) = E{F{x)f - E'^{F{x)) = E{F{x)f - F^{x) 

n — 1 


= 1 - 2 — 


^ J=0 


.kJl 


n-l-j 

X a'^r(n-j) ^ ■ 

i=0 


aln(l)] 


i-i± 


Finally 


, Q, n—1 

MSEiFix)) = -^^C(2n-2,jKr(n-j) [-ln(l) 


i=o 

[In (I)]’ 


E 

2=0 


2a 


(13) 
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2.1 The rth estimate of f{x) and F{x) 

To find the the rth estimate of f{x), we have 


Eif{x)Y = 


if{x)ygiw)dw = 


roo ^r^rw 


(any 


^riw+i) r(n)ii;”+i 




fc\™ 


,,n—r+l 


(an)” 

r(n)x’' 

(an)” 

^(n)x-^ Jo 

(any 


,,n—r+l 


dw 


■dw 


r°° “ 


E 


r(n)x” Jo ^ 


(In (-))'^ , 


(an)” .^rJ(ln(|)) 

r(n)'i’’' ^ 


j /*cx) ^—OLnjw 






,,n-r-j + l 


dw 


(«n)” y>r^'(ln(|)) 

r(n)'i’’' ^ 


J 




j=o 


J'- 


{any - r^(ln(|)) 

r(n)'i’’' 


t=o 




r(n - r - j) — 




n— 1 —1 

1 r-^ 


r(n)x’' 


E 

i=o 


(In(l)V 


J'- 


r(n — r — j)(an)^''"’’. 


Also, the rth estimate of F(x) can be found by calculating the following integral. 

(an)” 


E(F(x)y = 


/ (F(x)yg(w)dw = / 
'0 Jo 


1 - 




'0 j=0 


k 

X 

an) 


T(n)i 


,n+l 




T(n)i 


,n+l 


e-^^-^l^dw 


r(n) ^ Jo 


oo / l^\jw ^ — an/'u 


nn+1 


-dw 


[an 


r „oo jwln(^) -an/w 

rMyy’-X(-0’l - 


w 


■n+l 


(an)"" 

(an)” 

W) 






t=o 


2=0 

OO 


^C(r,j)(-1)^E 


d Jo 


n+l—2 


dw 


n— 2 —1 ^—anz 


dz 


t=o 


1=0 

OO 


^ ’ J=0 2=0 


(dln(f))’ 


Then 


n—1 






r(n) ^ 

^ ' 3=0 


T(n-i)(—y 

an 


r(n — f)(an)*. 


i=0 


(14) 


(15) 
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2.2 The rth estimate of f{x) and F{x) 
The rth estimate of f{x) is easily obtained as follows. 
r . 


E{f{x)r = {f{x)r 


(n — 1)! 


t-“-^[ln(t)-nln(/c)]"-Mf 




(n - lya^k 


an /*oo 


[ln(t) - ln(a;) - (n - 1) ln(fc)]'-(" 


x''{n — 1)! Jxk”—'^ [ln(t) — nln(fc)](’’“i)(”“i) 


Let z = ln(t) — nln(fc) ^ dz = jdt, then 


E{f{x)y = 


[z-ln(a;)+ln(fc)]’'("-2)A:-“"e- 


’(n-1)! 7ln(|) 


^(r-l)(n-l) 


-dz 


(n-iya^ r°° 
xy-n-iy. 

r T\r n rin-2) 

(n-iya^ y. ^ 
xyn-iy. ^ 

/ 1 \r n »’("-2) 

[n — 1) a 
x''{n — 1)! 


C(r(n - 2), j) (- In (f))' zriu- 2 )-j 




dz 


C(r(n-2),j) ( 


j=o 



z^-^-i-y-»^dz 


^ r(n-r-j) " ^ ' exp (-gin (I)) (aln(f))' 


a^-r-j 


t=0 


ll 


Therefore 


Eif{x)) = 3fT77::-TT7 C(r{n - 2), j)a^ 


x“+’'(n- 1)! 


i=o 


n-r-j-l i 


— In 


©]T(n-.-;) X: £±iM. , 16 ) 


2=0 


Also, the rth estimate of F(x) is similarly obtained as follows. 

EiF{x)y = [ {F{x)yh*{t)dt 


1 - 


[ln(t) — ln(x) — (n — 1) In(fc)]^ 


' xk'^~^ L 

/*00 E 


[ln(t) — nln(fc)]^ 


pxk^ 


h*it)dt+ Vxh*it)dt 


Ik" 


Ixk"-^ j=0 


Yc{r,j){-iy 


[ln(t) — ln(x) — (n — 1) ln(fc)]^(” 


[ln(t) — n ln(A:)]-l(”“i) 


pxk^ 


(n- 1)! 

^,nj^na. 


t “ ^[ln(t) — nln(fc)]" ^dt + 


Ik" (n-1)! 


- n\n{k)y-^dt 


a 


[ln(t) - ln(x) - (n - 1) ln(fc)]^(” 


[n - 




^y ^C{r,jy [ln(t)-nln(fc)]0-i)(n-i) 

nxk^~^ 

- t-°‘-y\n{t)-n\n{k)y-^dt. 

■ Jk" 


(n-1)! Jkn 
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Let z = ln(t) — n\n{k) ^ dz = \dt, then 


EiF^r = i;;T:7iTEC(''.J)(-l) 


[z — ln(x) + ln(fc)]'^(” 


Then 


+ 


+ 


3 R code 


j=0 


I Inn 




rMf) 


^n—1 i^ — an^—ocz 


dz 


1 )! Jq 

a" .. r C{j{n-l)n){-Hl)yz^ 

/ln(f) 






in — 1)! a" \k) ^ 

' ' L i=0 

r i("-i, 

]^C{r,j){-iy COXn-l),z) (-ln(£ 


-dz 


(n- 1)! 




i=o 


i=0 


/ln(f) 


1 - 


/xx-“!^'(aln(|))* 
\k) ^ i\ 


i=0 


i("-i) 


^ j = 0 2=0 


x\Y r(n — i) 


n—i—l 

E 

z=o 


2=0 

exp(-aln(f))(aln(f))' /x\-“ 


a'‘ 


l\ 


(f) (f) -E 


o -i(aln(|)r 


2=0 




i?(F(x))’- = 


(n — l)!x^ 


E C'0(^-0.0 (-ln(l)) aT(n-i) 


i=o 


2=0 


n— 2 —1 

E 

i=0 


(«ln(f)y , /fc 


n 


+ I “ 


« :^'(aln(|)r 


(I) -s 


2=0 


(17) 


The R code to compare the bias and MSE of the estimators is as follows. 

sim=function(t,n,k,alpha,r) 

{ 

sfh<-0 

sFh<-0 

sft<-0 

sFt<-0 

ford in l:t){ 
x<-array(, c(l,n)) 
for (i in l:n) { 

u<-runif(1,0,1) 

X[i] <-k*(1-u)~(-1/alpha)} 
alphah<-n/sum(log(x)-log(k)) 
fx<-alpha*k~alpha/x[1]~(alpha+1) 

intB<-function(z){ z"(-n)*exp(-alpha*n/z)*(k/x[1])~z} 
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B<-integrate(intB,lower=0,upper=Inf)$value 
Eftildx<-(alpha*n)~n/factorial(n-1)/x[1]*B 

intBl<-function(z){ z"(-n+1)*exp(-alpha*n/z)*(k/x[1])~(2*z)} 
Bl<-integrate(intBl,lower=0,upper=Inf)$value 
Eftildxs2<-(alpha*n)~n/factorial(n-1)/x[1] ~2*B1 
MSEftildx<-Eftildxs2-2*fx*Eftildx+fx~2 

intA<-function(z ){ z~ (2*n-4)*exp(-alpha*z)/(z+log(x[1])-log(k))"(n-1)} 
A<-integrate(intA,lower=0,upper=Inf)$value 

MSEfhx<-(n-1)*alpha"n*k~alpha/x[1]"(alpha+2)/factorial(n-2)*A-alpha~2* 
k"(2*alpha)/x[1]~(2*alpha+2) 

Fx<-1-(k/x[1])"alpha 

intB2<-function(w){ w"(-n-1)*exp(w*log(k/x[1]))*exp(-alpha*n/w)} 

B2<-integrate(intB2,lower=0,upper=lnf)lvalue 
EFtildx<-l-(alpha*n)"n/factorial(n-1)*B2 

intB3<-function(w){ w"(-n-1)*exp(2*w*log(k/x[1]))*exp(-alpha*n/w)} 
B3<-integrate(intB3,lower=0,upper=lnf)lvalue 

EFtildxs2<-l-2*(alpha*n)"n/factorial(n-l)*B2+(alpha*n)"n/factorial(n-l)*B3 
MSEFtildx<-EFtildxs2-2*Fx*EFtildx+Fx"2 

intAl<-function(w){ (w-log(x[1]/k))"(2*n-2)*exp(-alpha*w)*w"(-n+1)} 

Al<-integrate(intAl,lower=log(x[1]/k),upper=lnf)lvalue 

MSEFhx<-l-2*(k/x[1])"alpha+alpha"n/factorial(n-1)*A1-(l-(k/x[1])"alpha)"2 

sfh<-sfh+MSEfhx 

sft<-sft+MSEftildx 

sFh<-sFh+MSEFhx 

sFt<-sFt+MSEFtildx 

} 

mMSEfhx<-sfh/t 

mMSEftildx<-sft/t 

mMSEFhx<-sFh/t 

mMSEFtildx<-sFt/t 

return(c(mMSEfhx,mMSEf tildx,mMSEFhx,mMSEFtildx)) 

} 

simdO, 5,1; 5,1) 

siml=function(t,k.alpha,r){ 
i<-seq(3,35,1) 
for (j in i){ 
sim(t,j,k,alpha,r)}} 
siml(10,1,5,1) 


4 Tables 

In order to get the idea of efficiency between the two type of estimation i.e MLE and UMVUE. We have 
generated a sample of size 4(1)15(5)100 from the Pareto distribution with q;= 0.5(0.5)2 and A:=0.5(0.5)2. 
We have given Tables based on one thousand independent replication of each experiments. 

Table 1. shows the bias and MSE of the estimators of the pdf and bias and MSE of the estimators of cdf 
are shown in Tables 2. The value in the bracket is for the MSE in each tables. From the Tables, it has 
been seen that MLE of pdf and cdf are more efficient than UMVUEs. 

One should note that UMVUE of a is better than MLE of a. 
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Table 1. MSB of f{x) and f{x) for different values of a and k respect to n 


n 

Q = 0.5 
k = 0.5 

0 = 1 

k = 1 

Q = 1.5 
k = 1.5 

a = 2 

k = 2 

Q = 0.5 

k = 2 

a = 2 
k = 0.5 

4 

.4723690000 

(.4551780000) 

,4755806957 

(.4572692552) 

.476511 

(.457605) 

.476705 

(.457264) 

.029389 

(.028230) 

7,658800 

(7.374350) 

5 

,3172310000 

(.2725650000) 

,3207780187 

(.2757554306) 

,322473 

(.277487) 

.321924 

(.276551) 

.019827 

(.017036) 

5,157680 

(4.434020) 

6 

.2383780000 

(.1881800000) 

.2404067792 

(.1896539273) 

,241698 

(.190811) 

.242471 

(.191536) 

.014881 

(.011743) 

3,885290 

(3.071030) 

7 

.1914040000 

(.1421410000) 

,1929655509 

(.1432299232) 

,193791 

(.143884) 

.194454 

(.144461) 

.011901 

(.008828) 

3,105110 

(2.305270) 

8 

.1599510000 

(.1133700000) 

,1602595912 

(.1133783291) 

,161806 

(.114652) 

.161604 

(.114415) 

.009946 

(.007042) 

2,588020 

(1.832800) 

9 

.1372420000 

(.0937400000) 

.1380770369 

(.0942468652) 

.138780 

(.094764) 

.139051 

(.094956) 

.008601 

(.005877) 

2.213210 

(1.509400) 

10 

.1195130000 

(.0791650000) 

.1209664617 

(.0801522609) 

,121456 

(.080488) 

.121541 

(.080528) 

.007514 

(.004982) 

1.944830 

(1.288590) 

11 

.1063690000 

(.0687560000) 

.1075933345 

(.0695600949) 

.107972 

(.069808) 

.108015 

(.069818) 

.006661 

(.004307) 

1.728000 

(1.116900) 

12 

.0957350000 

(.0606160000) 

,0966263202 

(.0611705833) 

,096937 

(.061363) 

.097635 

(.061865) 

.005985 

(.003789) 

1.561640 

(.989450) 

13 

.0873630000 

(.0543730000) 

,0876848730 

(.0545292384) 

.088438 

(.055038) 

.088472 

(.055049) 

.005442 

(.003386) 

1.418000 

(.882550) 

14 

.0802730000 

(.0492200000) 

.0808533231 

(.0495643965) 

.081013 

(.049654) 

.081069 

(.049681) 

.005006 

(.003069) 

1.299260 

(.796410) 

15 

.0735700000 

(.0444890000) 

.0746156473 

(.0451395324) 

.074685 

(.045168) 

.074919 

(.045318) 

.004613 

(.002790) 

1.197330 

(.724140) 

20 

,0532991000 

(.0307735000) 

,0538278427 

(.0310780487) 

,053924 

(.031129) 

.054039 

(.031196) 

.003330 

(.001923) 

,865418 

(.499644) 

25 

.0418259000 

(.0234799000) 

,0420976064 

(.0236273808) 

,042261 

(.023720) 

.042273 

(.023724) 

.002612 

(.001466) 

.676949 

(.379933) 

30 

.0341756000 

(.0188200000) 

.0346115111 

(.0190619643) 

.034740 

(.019133) 

.034775 

(.019152) 

.002141 

(.001179) 

.556918 

(.306728) 

35 

.0290451000 

(.0157795000) 

.0294767054 

(.0160172951) 

.029444 

(.015996) 

.029448 

(.015996) 

.001821 

(.000989) 

.473159 

(.257090) 

40 

,0252562000 

(.0135817000) 

,0254701187 

(.0136957411) 

,025650 

(.013795) 

.025593 

(.013761) 

.001581 

(.000850) 

.411040 

(.221060) 

45 

.0224046000 

(.0119535000) 

.0225300956 

(.0120188221) 

.022630 

(.012073) 

.022658 

(.012088) 

.001395 

(.000744) 

.362473 

(.193366) 

50 

.0198892000 

(.0105411000) 

.0202331230 

(.0107253823) 

.020334 

(.010780) 

.020295 

(.010758) 

.001253 

(.000664) 

.324570 

(.172036) 

55 

,0181373000 

(.0095638000) 

,0183714478 

(.0096881397) 

,018323 

(.009661) 

.018382 

(.009692) 

.001130 

(.000596) 

,294111 

(.155077) 

60 

,0165965000 

(.0087132000) 

,0167227711 

(.0087790315) 

,016761 

(.008799) 

.016767 

(.008802) 

.001036 

(.000544) 

,269022 

(.141232) 

65 

.0152848000 

(.0079948000) 

,0153716098 

(.0080394109) 

,015452 

(.008082) 

.015489 

(.008102) 

.000954 

(.000499) 

,247479 

(.129436) 

70 

.0141436000 

(.0073741000) 

.0142593414 

(.0074341631) 

.014292 

(.007451) 

.014337 

(.007474) 

.000880 

(.000459) 

.228735 

(.119243) 

75 

.0131255000 

(.0068239000) 

,0132738528 

(.0069011346) 

,013345 

(.006938) 

.013349 

(.006940) 

.000820 

(.000427) 

,213970 

(.111250) 

80 

.0123438000 

(.0064025000) 

.0124114845 

(.0064370354) 

.012482 

(.006474) 

.012486 

(.006476) 

.000766 

(.000397) 

.199643 

(.103542) 

85 

.0115424000 

(.0059735000) 

,0117098555 

(.0060606327) 

.011725 

(.006068) 

.011738 

(.006075) 

.000722 

(.000374) 

.188222 

(.097417) 

90 

.0109372000 

(.0056499000) 

.0110042351 

(.0056841841) 

.011088 

(.005728) 

.011080 

(.005724) 

.000680 

(.000351) 

.177858 

(.091881) 

95 

,0102967000 

(.0053096000) 

.0104005835 

(.0053631148) 

,010482 

(.005405) 

.010493 

(.005411) 

.000644 

(.000332) 

,167683 

(.086471) 

100 

.0098069000 

(.0050495000) 

.0098648715 

(.0050790332) 

.009939 

(.005117) 

.009971 

(.005134) 

.000615 

(.000316) 

.159213 

(.081978) 


The figures in the bracket refers to the MSE of MLE of f{x) and without bracket refers to the MSE of UMVUE of f{x) {f{x)) 
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Table 2. MSE of F{x) and F{x) for different values of a and k respect to n 


n 

Q = 0.5 
k = 0.5 

0 = 1 

fc = 1 

Q = 1.5 
k = 1.5 

a = 2 

fc = 2 

Q = 0.5 

fc = 2 

0 = 2 
fc = 0.5 

4 

,1333908446 

(.0014763889) 

,1368410441 

(.0012562228) 

.038435 

(.001100) 

.115906 

(.000794) 

,121205 

(.001880) 

.157459 

(.002199) 

5 

,1118961062 

(.0072949912) 

,1799791911 

(.0356532489) 

,172372 

(.029683) 

.198062 

(.055201) 

,137497 

(.013072) 

,124812 

(.009786) 

6 

.1408169615 

(.0037596610) 

.1737093067 

(.0032023158) 

.101804 

(.001684) 

.188658 

(.005213) 

,143735 

(.003836) 

.134947 

(.005169) 

7 

,2191326212 

(.0229514291) 

,2065645792 

(.0165453192) 

,195758 

(.013057) 

.250334 

(.067911) 

,126943 

(.004378) 

,160016 

(.007159) 

8 

,2007151499 

(.0034259269) 

,2012409768 

(.0032392638) 

.195431 

(.004168) 

.219443 

(.000621) 

,197351 

(.002333) 

.203045 

(.007496) 

9 

.2576817407 

(.0275929634) 

.2512987261 

(.0205299145) 

.147364 

(.004077) 

.135785 

(.003471) 

.124131 

(.002910) 

.269592 

(.060899) 

10 

.1492978531 

(.0023002929) 

.1837933180 

(.0030831665) 

.150179 

(.003979) 

.220440 

(.005504) 

,134771 

(.001508) 

.213361 

(.007332) 

11 

.2623413820 

(.0139901431) 

.2259891179 

(.0076116749) 

.259195 

(.012447) 

.274008 

(.046732) 

.163600 

(.003893) 

.129209 

(.002419) 

12 

,1819439058 

(.0062366813) 

,1679376346 

(.0049937937) 

,197488 

(.004988) 

.176699 

(.004595) 

,147436 

(.004263) 

.205010 

(.005586) 

13 

.2177045578 

(.0058142184) 

.1674005776 

(.0033358805) 

.235847 

(.006945) 

.254445 

(.008379) 

,141234 

(.002349) 

.141272 

(.002350) 

14 

.2356903001 

(.0060898969) 

.0548103487 

(.0003127263) 

.227515 

(.006204) 

.193420 

(.004149) 

,188948 

(.006286) 

.136300 

(.001994) 

15 

.1076037454 

(.0011259936) 

.1555011068 

(.0024132218) 

.122969 

(.001482) 

.253257 

(.007139) 

.199792 

(.004118) 

.191199 

(.003743) 

20 

,2627256864 

(.0065292099) 

,2622509943 

(.0062066109) 

,248371 

(.005120) 

.232384 

(.004292) 

,164861 

(.002324) 

,168388 

(.002052) 

25 

.2488455377 

(.0041453670) 

.2598193829 

(.0052022850) 

,218770 

(.004877) 

.259375 

(.004883) 

,151849 

(.001275) 

.177311 

(.001795) 

30 

.2251732413 

(.0042208658) 

.2172670151 

(.0023800197) 

.257171 

(.004369) 

.211684 

(.003983) 

.113000 

(.000552) 

.118642 

(.000612) 

35 

.2458125666 

(.0028694984) 

.2072001232 

(.0018011252) 

.198559 

(.003254) 

.204366 

(.003352) 

.178507 

(.001266) 

.156018 

(.000935) 

40 

,1737078950 

(.0010316986) 

,2554091293 

(.0032353869) 

,251794 

(.002759) 

.152611 

(.002133) 

,116046 

(.001487) 

.156706 

(.000818) 

45 

.1940551896 

(.0011808928) 

.2429150460 

(.0021535385) 

.173837 

(.002250) 

.222621 

(.002867) 

.157868 

(.000733) 

.108663 

(.000325) 

50 

.1064113870 

(.0011050152) 

.2267825379 

(.0015767910) 

.179884 

(.000884) 

.228218 

(.002635) 

.125222 

(.000394) 

.151083 

(.000594) 

55 

,2299332896 

(.0024125373) 

,1855659085 

(.0008607060) 

,252088 

(.002056) 

.218652 

(.002338) 

,163238 

(.000639) 

,181364 

(.001967) 

60 

,1392457773 

(.0013392528) 

,1773309892 

(.0017724707) 

,127104 

(.001195) 

.227032 

(.002202) 

,170991 

(.000649) 

,108474 

(.000238) 

65 

.1986728058 

(.0018421105) 

.2322009300 

(.0012944937) 

.253793 

(.001897) 

.115857 

(.000989) 

,107674 

(.000215) 

.118045 

(.000262) 

70 

.2288138955 

(.0019021751) 

.1496013496 

(.0012705637) 

.251293 

(.001864) 

.108626 

(.000203) 

.134125 

(.000320) 

.123655 

(.000268) 

75 

,2524003249 

(.0017059572) 

,1656385542 

(.0013387601) 

,252058 

(.001537) 

.226321 

(.001769) 

,113991 

(.000209) 

,124311 

(.000252) 

80 

.2232959249 

(.0016486241) 

.1154621504 

(.0008181671) 

.238145 

(.001686) 

.123070 

(.000887) 

.159104 

(.000408) 

.104494 

(.000719) 

85 

,1569167081 

(.0003712149) 

.1116721695 

(.0007421171) 

.197303 

(.000637) 

.038435 

(.000164) 

,100752 

(.000141) 

.160423 

(.000959) 

90 

.2498099645 

(.0014551335) 

.1551047337 

(.0003408055) 

.154748 

(.000339) 

.084016 

(.000483) 

.080884 

(.000459) 

.071267 

(.000064) 

95 

,2525238629 

(.0013016097) 

,1423488455 

(.0002657662) 

,237976 

(.000956) 

.038663 

(.000151) 

,021378 

(.000005) 

,055160 

(.000255) 

100 

.1352354632 

(.0008108742) 

.0786840925 

(.0004021173) 

.105801 

(.000132) 

.064531 

(.000305) 

.055326 

(.000245) 

.026498 

(.000081) 


The figures in the bracket refers to the MSE of MLE of F{x) {F{x)) and without bracket refers to the MSE of UMVUE of F{x) {F{x)) 
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1. Introduction 


The Pareto distribution has been used in connection with studies of income, property values, 
insurance risk, migration, size of cities and firms, word frequencies, business mortality, service time 
in queuing systems, etc. 

The objective of this paper is to discuss efficient estimation of probability density function (pdf) 
and cumulative distribution function (CDF) of Pareto distribution which has been one of the most 
distinguished candidates for the honor of explaining the distribution of incomes, assets, etc. 

We assume that random variable X has Pareto distribution with parameters a and k (known) and 
its pdf is as 

ak“ 

f{x) — — 0 < k < X, a > 0, 
and CDF is 


Fix) = 1 - 


a 


k < X. 
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In economics, where this distribution is used as an income distribution, k is some minimum income 
with a known value. Asrabadi [1] derived the uniformly minimum variance unbiased estimator 
(UMVUE) of pdf, CDF and the rth moment. 

In general, we expect that the UMVU estimators are better than MLEs. In Pareto distribution, we 
show that UMVU estimators of parameter a is more efficient than the MLE, but for pdf and CDF, ML 
estimators are biased and more efficient than the UMVUEs. 

2. Maximum likelihood estimator 


Let Xi,X 2 ,... ,Xn be a random sample of size n from the Pareto distribution. According to ML 
method we obtain the MLE of a and it is given as a where a — n( ln(|))^'. 

Therefore, by using the property of MLE we can obtain the estimator of pdf and CDE with 
replacement of a instead of a in the pdf and CDE, respectively. Then 


fir 

/ (x) = , a > 0, 0 < k < X, 


F(x) = 1 - 



0 < k < X, fi > 0. 


( 1 ) 

( 2 ) 


We know that pdf of S = Yl'i=t 


g(s) = 


fin) 


exp(—as), s > 0, 


and by using some elementary algebra, we can find the distribution of w = fi as 


(an)'' 

g(w) = exp 

r(n)w''+‘ 


an i 

-[, w > 0. 

w 


Note. It Is clear that the MLE of a is biased and MS£(fi) = 
Theorem 1. (A) /(x) is a biased estimator off(x) and 

(any+^ 


(n-l)2(n-2)- 


■j 11 — z 

WW) = 7^ E 


J! 


r(n-j-l)(ln( 


(B) F(x) is a biased estimator of F(x) and 

I 

E(F(x)) = 1 


1 1^ (ar 


^“"^(n-i)(ln(f^ 


(3) 

(4) 

(5) 

( 6 ) 


Proof. In cases of (A), we can easily find the expectation of /(x) with substituting this formula: 
(^)iB _ gi«in(^) _ Also, the GAMMA function defines for variable grater than zero, 

then) must be less than (n — 1) and the proof is complete. In the case (B), the proof is similar as in the 
case (A). □ 

Theorem 2. 

(A) 


MS£(f(x)) = 


1 

— y 


I^anf^^ / fk 

. . r(n - J - 2) In - 

r(n)x^ j'. V 


ar 


r(n)x‘ 


•a+2 




j=0 


17 


\x‘ 


(7) 
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(B) 


MS£(f(x)) = 


1 


2^ {any ( (k 

rt - 

r{n) ]'■ \ V 




( 8 ) 


Proof. In cases (A) and (B) we should find E(f (x))^ and E(F(x))^ as the previous Theorem, respectively. 
So by using some elementary algebra the proof is complete. □ 

3. MSE of UMVU estimator 

Asrabadi [1] derived the UMVUE of a,f{x) and E(x). Here UMVUE of a,f(x) and F(x) are denoted 
by a,f{x) and F(x), respectively. So 
n - 1 


ln(t) - n ln(k) 

(n - l)[ln(t) - ln(x) - (n - 1) \n{k)f-^ 


and 




Fix) = 1 - 


x[ln(t) - n ln(/<)]"-i 
[ln(t) — ln(x) — (n — 1) ln(k)]"^’ 


[ln(t) — n ln(k)]" ' 
where k <x < tk^^", and t = x, is the observed value of T. 

Theorem 3. 

(A) 

(n-l)a2fe“ "-3 


MSE (fix)) = 


^ C(2n - 4J)aJr(n - j - 2) In 

j=o 


Fin - l)x“+2 


"-3-i 


E 

1=0 


(ln(j))' / ak 


( x“+i 


(B) 


MSF(F(x)) = 


k“ 

r(n)x“ ^ 


^ C(2n - 2J)a'rin -j) (^- In 


X i: 

1=0 


a'(ln(f))' 


(9) 

( 10 ) 

( 11 ) 


( 12 ) 


(13) 


whereCin,k)^^^. 

Proof. In cases (A) and (B). we can obtain E(f(x))^ and E(F(x))^ by using pdf of T that is given in [1]. 
In the process to calculate the integral we should note that 


f 


Fin) 


^ a 


Hence, the proof is complete. □ 

2 "I ft 

Note. One should note that MSF(a) = 
















MSE 
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Fig. 1. Comparison of the MSE of the estimators of pdf and CDF with respect to observation generated from the Pareto 
distribution. „ 
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4. Comparison of MLE and UMVUE 

It is obvious that the UMVU estimator of a is more efficient than the MLE for any value of n. Now 
in order to get the idea of efficiency between MLE and UMVUE of pdf and CDF, we have generated a 
sample of size 4(1)15(5)100 from the Pareto distribution with a = 0.5(0.5)2 and k — 0.5(0.5)2. We 
have given graphs based on one thousand independent replications of each experiments (Fig. 1). From 
the graphs, it has been seen that MLE of pdf and CDF are more efficient than UMVUEs. 
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results indicate that there are no absolute superiorities of MLEs 
over the UMVUEs of PDF and RF and vice versa. 
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1. Introduction 


We consider a random variable X has the Pareto distribution (PD) with PDF 




ak“ 
x“+’ ’ 


and RF 


( 1 ) 


G(x) = Prob{X > x} = 



( 2 ) 


where a is a shape parameter (a > 0), and k is a scale parameter (known, and x > k > 0). in 
economics, k usually represents some minimum income with a known value, see [1]. 

PD was applied by Pareto [7] to model the allocation of wealth among individuals and the dis¬ 
tribution of incomes. It has been widely used in economics, insurance (general liability, commercial 
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auto [9]), geography (sizes of human settlements [8]), physical sciences (sizes of sand particles or 
meteorites [8], clusters of Bose-Einstein condensate near absolute zero [5]), chemical sciences (dis¬ 
tributions of electrolytic powder production [4]). Asrabadi [1] established the UMVUEs for the PDF 
and cumulative distribution function (CDF) of PD. Based on the work of Asrabadi [1], Dixit andjabbari 
Nooghabi [2] tried to study the mean square errors (MSEs) of the MLEs and UMVUEs for the PDF and 
CDF of PD and their results seem to show that the MLEs are more efficient than the UMVUEs of PDF 
and CDF. Unfortunately, their work are seriously flawed. Most of their main claims in [2] are wrong, 
and their conclusion, the MLEs are more efficient than the UMVUEs of PDF and CDF, is unreasonable. 

We present our main results in Section 2. Most of the results in Section 2.1 are corrected versions of 
the wrong results of [2]. We also notice that the exact expressions of the MSEs of estimators of PDF and 
RF may not be useful in case of large scale samples and large scale numerical computations. For this 
reason we have derived the asymptotic expressions of the r-th raw moments and MSEs in Section 2.2. 
Two numerical examples are provided in Section 2.3 to show how to compute the efficient estimations 
of PDF and RF. In Section 2.4 we expose the fatal errors in [2]. 

2. Main results 

As a notational convenience, let 

X dz 1 

z = z(x) = log-, Z;, = Zx(x) = — = 
k ax X 

through the rest of this paper. It is known that the UMVUEs of/(x) and G(x) are given by [1] 
n — 1 / z\ "-2 „ , z\ "-1 

f(x) = ^\~7) ’ 

where z < s, s — z(Xi) and s follow the Gamma distribution Ga(n, a). Note that the UMVUE of 
O' is a = (n — l)/s. 

The MLEs of/(x) and G(x) can be computed easily, they are 
/(x) = aZxe~“, G(x) = 

where a = ^ is the MLE of a. Note that the PDF of s is given by h(s) = exp(—as). 

2. 1. The r-th raw moments of estimations 


Theorem T. For n > r > 0, the r-th raw moments off{x) and G(x) are given by 

E(J{x)y = (nazj’’ i^nraz)"^''Kn-r {2-Jnraz ), 
r(n) 

E(G(x)y = j^(^nrazyK„ {l^nraz ), 
where /^^(x) is the modified Bessel function [6]. 

Proof. For the proof we just need to note the well-known integral representation [6], 


/<v(x)= 




x^ \ dt 




(3) 

(4) 


Corollary 1. The mean square errors of f(x) and G(x) are given by 
MSE(/(x)) = (2V2naz^ 


r{n) 

4naZj 

'~r{n) 


f(x)(^/nazy (2Vnaz) -F/^(x). 


(5) 
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MSE(G(x)) = -j^(V2naz)'’Kn (ly/lnaz^ 


r(n) 


G(x)(y/naz)''I<„ (l^/naz) + G^(x). 


( 6 ) 


Theorem 2. For n > r > 1, the r-th raw moments of f(x) and C(x) are given by 


E (jf (x)/= (az,(n-I)/ Y(x) 




r(n-l) 


U(nr — n — r+1,r — n + 1, az), 


- r r(nr —r+1) 

£(G(x)) =-- -G(x)U(nr — n — r+1,1 — n, az). 


(7) 


( 8 ) 


where U(a, b, c) is the Kummer confluent hypergeometric function [6]. 

Proof. Note that the Kummer confluent hypergeometric function has an integral representation [6], 


U(a, b, c) = 


pC eoc 

— / t“-’ 

r(a) Jo 


(1 + t)'’ " ’ e "dt, 


r{a) JO 

and the proof is completed by applying the Kummer transformation [6], 
U(a, b, c) — c^^'’U(1 + a — b,2 — b, c). □ 


Corollary 2. The mean square errors of f (x) and G(x) are given by 

MSE(f(x)) = azj(x)(n- 1) ^G(n - 1,3 - n, az) -f(x). (9) 

r(n - 1) 

r(2n-l) , 

MSE(G(x)) = G(x)—- -U(n - 1, 1 - n, az) - G^(x). (10) 

r(n) 

Figs. 1 and 2 illustrate the efficient estimators between MLEs and UMVUEs of PDE and RE for 
n = {5, 6, 30}, /< = {1, 5), 0 < a < 10 and k < x < 50. In the graphs, the black areas indicate that 
the MLEs of PDE and RF are more efficient than the UMVUEs while the white areas mean the UMVUEs 
of PDF and RF are more efficient than the MLEs, and there is no evidence that the black areas or the 
white areas will disappear from the first quadrant. Thus we conclude that the MLEs are not generally 
more efficient than the UMVUEs of PDF/RF and vice versa. We also notice that Corollaries 1 and 2 can 
help us to obtain more efficient estimations, see Example 1. 


2.2. The convergence rate of estimators 

Corollaries 1 and 2 can be expediently used to calculate the MSEs of the estimations for a small 
sample. Flowever, in practice, we find that the corollaries would not be expedient for a large sample 
and the large-scale numerical computation. Eor reasonable large n, direct numerical evaluations of 
r (2n — 1), K„ (2y/naz) , U (n — 1,1 — n, az) will incur either overflow or underflow. Eor large-scale 
numerical computation, symbolic computations will run for a very long time. Therefore, it is necessary 
to study the large sample properties of the r-th raw moments and the asymptotic behaviors of the 
MSEs. 


Lemma 1. For fixed y >0, a > 0, then 


_ 00 

K„ {2^y(v + a)'^ = {v {2vy-)-^Y, 


^^ ai{a,y) 

1=0 


v 


( 11 ) 














52 


H. He et al. / Statistical Methodology 21 (2014) 49-58 





0 2 4 6 8 10 


a 

(e) Graph of MSE(/) - MSE{/) < 0 when 
n = 6 and k = 5. 



(f) Graph of MSE(/) - MSE(/) < 0 when 
n = 30 and k = 5. 


Fig. 1. Efficient estimators between MLE and UMVUE of PDF. 


as p ^ + 00 , where the first three coefficients are given by 

ao(a,y) = 1, ai(a,y) =y^/2 - (a+ l)y + (6a^ + 1)/12, 
a 2 (a,y) = y'‘/8 - (3a + 7)yV6 + 48a + 6l)yV24 

- (I44a^ + 144a^ + 312a+312)y/288 + (36a'‘ - 96a^ + Ua^ + l) /288. 
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0 2 4 6 8 10 


(a) Graph of MSE(C) - MSE(C) < 0 when (b) Graph of MSE(G) - MSE(G) < 0 when 
n = 5 and ff = 1. n = 6 and k= 1. 




(c) Graph of MSE(G) - MSE(G) < 0 when (d) Graph of MSE(G) - MSE(G) < 0 when 
n = 30 and ff = 1. n = 5 and k= 5. 




02468 10 02468 10 


(e) Graph of MSE(G) - MSE(G) < 0 when (f) Graph of MSE(G) - MSE(G) < 0 when 


n = 6 and k = 5. 


n = 30 and k= 5. 


Fig. 2. Efficient estimators between MLE and UMVUE of RE. 


Lemma 2. For fixed y >0, r>l, +oo>a> —oo, then 

r(nr — n — r + l)U(nr — n — r + 1, ar — r — n + 1,y) 

= (r - bi(a,r,y) ^ 


( 12 ) 


asn + 00 , where the first coefficient is given as bn(a, r,y) = 1, and when r — 2, the second and third 
coefficients are given by 

b, (a, 2, y) = + 2(1 - a)y + (Sa^ - 20a + 17) /8, 
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b 2 (a, 2,y) = yV2 - 2ay^ + (384a^ - 192a + 144)yVl28 
+ (-256a^ + 384a^ - 416a + 288) y/128 
+ (64a'‘ - 192a^ + 320a^ - 408a + 289) /128. 

Note that the proofs of Lemmas 1 and 2 are very lengthy, so we only provide the results here. For 
interested readers, please see the related sections on Laplace’s method in [10], 

Theorem 3. If n oo, then the following formulas hold: 

(a) £(f(x)y =/'-(x) + 0(n-’). 

(b) £(G(x))'- = C''(x) + 0(n-i). 

(c) E(fix)y =/''(x) + 0(n^’). 

(d) £(G(x))'' = C''(x) + (9(n-i), 

Proof. Cases (a) and (b) are obtained by applying Lemma 1 and Stirling's formula [6], the proofs for the 
cases (c) and (d)are similar to the cases (a) and (b) except applying Lemma 2 instead of Lemma 1. □ 

Theorem 4. For n^ oo, we have the following results: 

(a) MSE(f(x)) =/^(x) ((az - + pyaz)n-^ + 0(n^^)), 

(b) lVlSE(G(x)) = G^(x) ((az)^n^^ + qi(az)n^^ + 0(n^^)), 

(c) MSE(f(x)) =p{x) {{az - ifn^^ + p 2 (az)n-^ + 0(n^^)), 

(d) MSE(G(x)) = G^(x) ((az)^n^^ + q 2 {az)n^^ + 0(n^^)), 

where the functions pi,p 2 , qt, q 2 are defined as 

p, (x) = 7 xV 4 - 14x^ + 32x^ - 24x + 5, q,(x) = 7x^4 - 7x^ + 5x^ 

P 2 (x) = x^/2 — 4x^ + lOx^ — 8x + 2, q 2 (x) = x^/2 — 2x^ + 2x^, x > 0. 

Proof. The proof of Theorem 4 is similar to the proof of Theorem 3. □ 

Next, we will discuss the efficient estimations under the large scale sample. By Theorem 4, we have 
MSE(f(x)) - MSE(f(x)) ^ (f(x)/nfp(az), 
and 

MSE(G(x)) - MSE(G(x)) ^ {G{x)/nf q{az), 

where p{az) — pi(az) — P 2 iaz) and q(az) = qi(az) — q 2 (az). Note thatp(Q!z) = 5(az)"'/4 — 
lO(az)^ + 22(az)^ — 16(az) + 3, q(Q!z) = 5 (q!z)"'/ 4 — 5(az)^ + 3(az)^. We notice that the algebraic 
signs of p(az) and q(az) can approximatively determine the efficient estimators of PDF and RF when 
n is large. 

Corollary 3. For sufficiently large n, we have the following results: 

1. If az e [0.2856, 0.9168] U [1.857,4.9439], then p(az) < 0 and the MLF is more efficient than the 
UMVUF of PDF. 

2. If az e (0, 0.2856) U (0.9168,1.8537) U (4.9439, oo), then p{az) > 0 and the UMVUF is more 
efficient than the MLF of PDF. 

3. If az e [0.7351, 3.2649], then q(az) < 0 and the MLF is more efficient than the UMVUF ofRF. 

4. If az e (0, 0.7351) U (3.2649, oo), then q(c^ > 0 and the UMVUF is more efficient than the MLF of 
RF. 
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2.3. Numerical example 

Example 1. Efficient estimator in the small sample. 

We use Dyer [3] annual wage data (In multiples of 10,000 US dollars) to illustrate our results. The 
values of the data are given below: 

1.12, 1.54, 1.19, 1.08, 1.12, 1.56, 1.23, 1.03, 1.15, 1.07,1.25, 1.19, 1.28, 1.32, 1.07, 

1.51, 1.03, 1.04, 1.16, 1.40, 1.08, 1.05, 1.58, 1.04, 1.19,1.11, 1.01, 1.57, 1.12, 1.15. 

Here we suppose that the minimum wage is 10,000 US dollars. Then the pertinent data are n = 30, 
k = 1 and O' 5.4025 (the UMVUE of a). 

Further, let 

MSE(/(x)) - MSE(f(x)) = 0, 
then, the positive roots are 

Xi = 1.0512, X 2 = 1.1794, X 3 = 1.3825, X 4 = 2.5563. 

Therefore, if x e [1.0512,1.1794] U [1.3825,2.5563], /(x) is more efficient than/(x); if x e 
(1,1.0512) U (1.1794,1.3825) U (2.5563, 00 ), /(x) is more efficient than/(x). Similarly, let 

MSE(G(x)) - MSE(C(x)) = 0, 

then, the positive roots are 

Xi = 1.1384, X 2 = 1.8443. 

Hence, when x e [1.1384,1.8443], G(x) is more efficient than G(x); when x e (1,1.1384) U 
(1.8443, 00 ), G(x) is more efficient than G(x). 

Example 2. Efficient estimator in the large sample. 

To compare with our first example, we let k—1 and a — 5 .4025 (the values of k and a are the same 
as in the first example) In Corollary 3 to see that, when X e [1.0543,1.1850] U [1.4093, 2.4970], /(x) 
is more efficient than/(x): when X e (1,1.0543) U (1.1850,1.4093) U (2.4970, 00 ), /(x) is more 
efficient than/(x); if x e [1.1458,1.8300], G(x) is more efficient than G(x): if x e (1,1.1458) U 
(1.8300, 00 ), G(x) is more efficient than G(x). 

2.4. Some comments on [2] 

The errors of the main results of [2] can be seen clearly from the following simple numerical 
calculation. Let n = [5,6,30], k = [1,5], a = [1,5], and x = (2, 4, 6 , 8 ], the mathematical 
expectation values and MSE values are listed in Table 1, and some of them are negative which is clearly 
absurd. It is not hard to see Theorems 1-3 of [2] are all wrong, where their Theorem 1 is about the 
mathematical expectation expressions of/(x), f(x); their Theorem 2 is about the MSE expressions 
of/(x), F(x). and their Theorem 3 is on the MSE expressions of/(x), F(x). Furthermore, all of our 
numerical simulations (Figs. 1. 2, Examples 1 and 2) also show that the main conclusion of [2], the 
MLEs are more efficient than the UMVUEs of PDF and CDF, is false. 

3. Conclusion 


We have studied the efficient estimation in PD in our work. Our results show that the efficient 
estimations of PDF and RF of PD depend on f®(r variables (n, k, a, x). Letgi(x),..., g„,(‘^) denote the 
different estimations (MLE, UMVUE, Bayesian estimation, etc.) of PDF or RF, we construct the following 
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Table 1 


Numerical values of £(/(x)), £(f(x)), MSE(/(x)), MSE(f(x)), MSE(f(x)) and MSE(f(x)), forn = {5, 6, 30}, fc=jl.5}, a = 
{1,5} andx = {2, 4, 6, 8}. Note that the values are calculated hy using the results of Dixit and Jabbari Nooghabi [2|. 


(n, k, a, x) 

E(I) 

£(£) 

MSE(/) 

MSE(f) 

MSE(/) 

MSE(f) 

(5.1,5.2) 

-3.88e+2 

-1.29e+2 

3.76e+3 

2.25e+3 

1.41e+l 

5.73e+0 

(6,1,5.2) 

8.28e+2 

2.23e+2 

-2.11e+4 

-7.75e+3 

-5.31e+l 

-1.43e+l 

(30.1,5.2) 

8.51e-2 

9.67e-l 

-2.01e+7 

-1.89e+5 

-2.37e+l 

-2.11e-l 

(5,1,5.4) 

-1.67e+3 

-2.26e+3 

3.80e+3 

3.79e+4 

4.63e-l 

3.05e+0 

(6.1,5,4) 

7.14e+3 

7.77e+3 

-4.34e+4 

-2.61e+5 

-3.50e+0 

-1.53e+l 

(30,1,5,4) 

9.91e+5 

1.89e+5 

-7.1e+14 

-l.le+14 

-2.62e+7 

-3.76e+6 

(5,5,1,6) 

1.55e-l 

1.97e-l 

2.94e-2 

1.65e-2 

5.45e-3 

7.17e-3 

(6,5,1,6) 

1.54e-l 

1.92e-l 

-4.54e-3 

5.45e-3 

2.23e-3 

5.14e-3 

(30,5,1,6) 

1.42e-l 

1.71e-l 

4.64e-4 

8.71e-4 

4.49e-4 

8.14e-4 

(5,5,1,8) 

4.94e-2 

3.95e-l 

1.34e-l 

5.05e-l 

9.06e-3 

4.61e-2 

(6,5,1,8) 

8.71e-2 

4.18e-l 

-9.79e-2 

-2.13e-l 

-4.11e-3 

9.06e-3 

(30,5,1,8) 

7.85e-2 

3.83e-l 

5.29e-5 

3.07e-3 

5.76e-5 

2.99e-3 


estimation 

MSE(gi(x)) < MSEte(x)) 

MSE(g„,(x)) < MSE(g,(x)) 

where the parameters n,k,a are given. As an estimator of PDE or RE. it is more efficient than all of 
gi(x), i— 1 ,..., m. It is also clear thatg(x) may have discontinuities. 
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Appendix. Proofs associated with Theorems 1-4 and Lemmas 1-2 
Proof of Theorem 1. 


gW = 


g,(x), ifxe|x 


gm(x), ifxejx 


n 

1=1 

m 

n 

1=1 


pOO ^ pOO J. r 

E(f(x)y = j (f(x)yh(s)ds = J 

I (asy-’’-'' exp (-as - — ) d(as) 
f 

/ t exp —t - 

Jo V 4t 


exp(—as) 

r{n) 


ds 


(anz,y 
r(n) 
(cynz,y 
r{n) 

2 


r(n) 


(nazy (^nrazy {2^/nmz). 


Eq. (4) can be proved similarly, and the proof of Theorem 1 is complete. 

Proof of Theorem 2. 


E(i{x)y 


(f(x)yh(s)ds = (l - J)" aV- 


exp(—as) 

r{n) 


ds 


(a(n - -l)z,y 
r(n) 



; aZ\(n-2)r 

1-1 exp(—as)d(as) 

as 2 
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(a(n - \)Zxy 
r(n) 

iajn - l)z,)^ 
r{n)^^ 
(a{n - -l)z,y 
r{n)^^ 
(a(n - l)zj'' 


pOO 

/ _ Q;z)("-2)r exp(-t)dt 

Jaz 
poo 

/ it + azf-W-r)^(n-2)r exp(-t)dt 

Jo 


r 

iazr-' / 

Jo 


(1 + f) 


(n-l)(l-r).(n-2)r 


exp(—azt)dt 


r(n)e‘ 

= (az,(n- 

Eq. (8) can be proved similarly. 


(az)" ''r(nr -2r + 1)U(nr - 2r + 1, n - r + 1, az) 


,_,,^^r(nr-2r+l) 


r(n-l) 


U(nr — n — r + 1, r — n + 1, az). 


Proof of Lemma 1. 

(2vy^) = f exp (-, - >^) d, 

(y(v + a))2 /■“ / , y°\ / 

= -^-i, exp -j exp (-W (log t+-))-, 

where y > 0, a > 0 and v > 0. Let/(t) = log t + y/t, then 
/'(t) = t-y/t^ /"(t) =-t - 2y/t^ 

it is clear/(t) has a unique minimum 1 + logy at t = y. Then, we obtain Eq. (11) by Laplace’s 
method [10]. Then the proof of Lemma 1 is complete. 

Proof of Lemma 2. 


r{nr — n — r + l)U(m — n — r + l,ar — r — n + 1,y) 

fnr-n-r(i ^ ^^ar-nr-1 ^-yt 
(1 -L 

f 

_ f ^-n((l-r)log(t)+rlog(l+t)) 

Jo 

where y >0, r>l, +oo>a> —oo. Let/(t) = (1 — r) log(t) + r log(l + t). then 
/'(t) = r(l + t)-'-(r-l)r\ 

It is clear/(t) has a unique minimum rlogr — (r — 1) log(r — 1) at t = r — 1. Then, we obtain 
Eq. (12) by Laplace’s method [10]. Then the proof of Lemma 2 is complete. 

Proof of Theorem 3. In the cases (a) and (b), the results can be easily obtained by applying Lemma 1 
and Stirling’s formula 



OO 

r(n) = 

k=0 


Q 

n'’ 


Co = 1 , 



C2 = 


1 


In the cases (c) and (d). the results can be derived by using Lemma 2 and Stirling’s formula. 

29 

Proof of Theorem 4. The proof process of Theorem 4 is similar to the proof of Theorem 3. 
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